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Abstract: The fuel-optimal control problem arising in noncoplanar orbital 
transfer employing aeroassist technology is addressed. The mission involves 
the transfer from high Earth orbit to low Earth orbit with plane change. The 
complete maneuver consists of a deorbit impulse to inject a vehicle from a 
circular orbit to elliptic orbit for the atmospheric entry, a boost impulse 
at the exit from the atmosphere for the vehicle to attain a desired orbital 
altitude and finally a reorbit impulse to circularize the path of the 
vehicle. In order to minimize the total fuel consumption, a performance index 
is chosen as the sum of the deorbit, boost, and reorbit impulses. Application 
of Pontryagin minimum principle leads us to a nonlinear, two-point, boundary 
value problem, which Is solved by using a multiple shooting method. 

Nomenclature 

A ■ C Sp H /2m 

1 DO r s a' 

A = C Sp H /2m 

2 LR s a' 

b = R /H 

a • a 

C D : drag coefficient 

C : zero-lift drag coefficient 

DO 

C L : lift coefficient 

C : lift coefficient for maximum lift-to-drag ratio 

LR 

D : drag force 

E : maximum value of L/D 

m 

g : gravitational acceleration 

H : altitude 

H : Hamiltonian 

i : inclination 

J : performance index 
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K : induced drag factor 

L : lift force 

m : vehicle mass 

R : distance from Earth center to vehicle center of gravity 

: radius of atmospheric boundary 

R : radius of low Earth orbit 

c 

R : radius of high Earth orbit 

d 

R : radius of Earth 

E 

S : aerodynamic reference area 

t : time 

V : velocity 

v : normalized velocity 

£ : inverse atmospheric scale height 

7 : flight path angle 

: heading angle 

<r : bank angle 

0 : down range angle 

<f> : cross range angle 

5 : normalized density 

\ : costate (Langrange) variable 

jLt : gravitational constant of Earth 

* " C L /C LR 

p : density 

t : normalized time 

AV : characteristic velocity 

Av : normalized characteristic velocity 
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Subscripts 


c : circularization or reorbit at LEO 

d : deorbit at HEO 

e : entry to atmosphere 

f : exit from atmosphere 

s : surface level 

I. INTRODUCTION 

The main function of the space transportation system is to deliver 
payloads from Earth to various locations in space. Until now, this function 
has been performed by various rockets, the space shuttle, and expendable 
upper stages using solid or liquid propellants. In particular, considering 
the economic benefits and reusability, an orbital transfer vehicle (OTV) is 
proposed for transporting payloads between low Earth orbit (LEO) and high 
Earth orbit (HEO). The two basic operating modes contemplated for OTV are a 
ground-based OTV which returns to Earth after each mission and a space-based 
OTV which operates out of an orbiting hanger located at the proposed Space 
Station. 

In a typical mission, a space-based OTV, which is initially at the space 
station orbit (SSO), is required to transfer a payload to geosynchronous 
Earth orbit (GEO), pick up smother payload, say a faulty satellite, and 
return to rendezvous with the orbiting hanger at SSO for refurbishment and 
redeployment of the payload. The OTV on its return journey from GEO to SSO 
needs to dissipate some of its orbital energy. This can be accomplished by 
using an entirely propulsive (Hohmann) transfer in space only or a 
combination of propulsive transfer in space and aeroassisted maneuver in the 
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atmosphere. It has been established that a significant fuel savings and hence 
increased payload capabilities cam be achieved with propulsive and 
aeroassisted mauieuvers instead of all -propulsive maneuvers 1 . This leads to am 
aeroassisted orbital transfer vehicle (AOTV), which on its return leg of the 
mission, dips into the Earth’s atmosphere, utilizes atmospheric drag to 
reduce the orbital velocity and employs lift and bank angle modulations to 
achieve a desired orbital inclination. Basically, the AOTV performs a 
synergetic maneuver, employing a hybrid combination of propulsive maneuver in 
space and aerodynamic maneuver in the atmosphere. 

It is believed that the concept of aeroassisted orbital transfer opens 

new mission opportunities for the space transportation system, especially 

with regard to the establishment of the permanent space station. Fig. 1 shows 

the space transportation architecture relevant to aeroassist technology. The 

2-7 

optimization of fuel is an important aspect of orbital transfer missions. 

In this paper, we address the fuel-optimal control problem arising in 
noncoplanar orbital transfer employing aeroassist technology. The maneuver 
involves the transfer from HEO to LEO with a plane change and at the same 
time minimization of the fuel consumption. It is known that the change in 
velocity, also called the characteristic velocity, is a convenient measure of 
fuel consumption. For the minimum-fuel maneuver, the objective is then to 
minimize the total characteristic velocity for deorbit, boost, and reorbit 
(or circularization) for a specified change in inclination angle. Application 
of Pontryagin minimum principle leads us to a nonlinear, two-point, boundary 

value problem (TPBVP), which is solved by using a multiple shooting 

... 8-10 
method. 
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II. BASIC EQUATIONS 


For the orbital transfer problem, the following assumptions are made. 

(1) The initial HEO and final LEO orbits are circular, (ii) The mission is 
comprised of three impulses, (iii) The vehicle is represented as a constant 
point mass during atmospheric pass, (iv) A Newtonian inverse square 
gravitational field is used, (v) Earth’s rotation is neglected, (vi) The 
atmosphere is exponential, (vil) The vehicle has a parabolic drag polar. 

The complete mission from HEO to LEO with atmospheric pass is depicted 
in Fig. 2. It is composed of three impulses: first, a deorbit impulse AV at 

d 

HEO to inject the vehicle into a HEO-entry elliptic orbit, second, a boost 

impulse AV at the exit from the atmosphere for the vehicle to attain 
b 

sufficient velocity to travel along an exit-LEO elliptic orbit, and finally, 
a circularizing impulse AV_ to circularize the path of the vehicle. 

Consider the basic equations of motion for different phases of deorblt, 
aeroassist (or atmospheric flight), boost and reorbit (or circularization). 

Deorbit 

Initially, we assume that the spacecraft is in a circular orbit of 
radius R , well outside the Earth’s atmosphere, moving with a circular 

d 

velocity V = Vyt/R . Deorbit is performed by means of an impulse AV , to 

d d d 

transfer the vehicle from the circular orbit to elliptic orbit with perigee 
low enough to intersect the dense part of the atmosphere [Fig. 2). Since the 
elliptic velocity at D is less than the circular velocity at D, the impulse 
AV is executed so as to oppose the circular velocity V . The deorbit impulse 

d d 

AV causes the vehicle to enter the atmosphere of radius R with a velocity 
d a 

V^ and flight path angle It is known that the optimal -energy loss 
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maneuver from the circular orbit is simply the Hohmann transfer and the 
impulse is parallel and opposite to the instantaneous velocity vector. 


Using the principle of conservation of energy and angular momentum at 
the deorbit point D, and the atmospheric entry point E, we get, 11 

V 2 /2 - M/R ■ (V-AV ) Z /2 - \>J R (1) 

e / a d d ' d 

RVcos(-y ) = R (V-AV ) (2) 

a e e d d d 

from which solving for AV . we get 

d 



2p ( 1 /R a - 1 /R d ) / [( R d /R a ) 2 /cos 2 r e - 1 ] 


(3) 


It is easily seen that the minimum value of the deorbit impulse AV 
obtained at y^ = 0, corresponds to an ideal transfer wherein the space 
vehicle grazes along the atmospheric boundary. To ensure proper atmospheric 
entry, the deorbit impulse AV must be higher than the minimum deorbit 

d 

impulse AV which is given by 

dm 

AV = T - J2p(l/R -1/R )/[(R /R ) 2 -ll (4) 

dm dH a d' L da J 

Aeroassist (Atmospheric) Flight 

During the aeroassist (or atmospheric) flight, the vehicle performs a 
three-dimensional skip maneuver and using aerodynamic lift and bank angle 
achieves the necessary the plane change. In this process, the vehicle 
decelerates due to the atmospheric drag. 

The equations of motion for the vehicle during the atmospheric pass are 

2 

given below (Fig. 3). The kinematic equations are, 
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§ - Vsln * 

^ = Vcosycos^/Rcos^ 
at 

^ = Vcosysin^/R 
at 

The force equations are 

dV _ , 

m-rr = - D - mgsiny 


(5a) 

(5b) 

(5c) 

(5d) 


mV^ = Lcos«r + m(V?R - g)cosy 

Qt 

“ Lsintr/cosy - (mV?R)cosycos0tan0 
dt 

where, 

L = C^pSV 2 /2; D = C D pSV 2 /2; C D - + KC 2 

g = p/R 2 ; R = H + R £ ; p = p^expt-H/S) 

Using the normalized variables, 

t = t /\ R^/p ; v = V/v'pTIT 

and the dimensionless constants, 

h = H/H ; b = R /H ; 5 = p/p = exp(-h£H ) 

a a' a s a 

" * SAb 1 C u,-^ 


in (5), we get the normalized form as 


dh 

dx 


= bvslny 


(5e) 
(5f ) 


( 6 ) 

(7a) 

(7b) 

(8a) 
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(8b) 


d0 _ bvcosycosft 
dx (b-l+h)cos^ 


d <f> _ bvcosysin 
dx " (b-l+h) 


(8c) 


£ - - A b(l+VW - b8sl " r 

dT 1 (b-l+h ) 2 


(8d) 


dy 

^ = A 2 bT)5vcos<r + 


bvcosy 

(b-l+h) 


,2 

b cosy 
(b-l+h)v 2 


(8e) 


di/i _ A 2 b6T > vsln<r ^ bvcosycos^itan# 
dx ~ cosy (b-l+h) 


where, A — CSp H /2m; A = C Sp H /2m 

1 DO s a' 2 LR s a' 

From the above equations of motion, we see clearly that during the 

atmospheric maneuver, if the lift vector L is rotated about the velocity 

vector V through the bank angle <r, it creates a lateral force component Lsin<r 

orthogonal to the vertical plane that has the effect of changing the heading 

angle ifi. At the end of the maneuver, the vehicle is already in vacuum and 

hence there is no lift. The equations (5c) and (5f) for the cross range angle 

2 

<p, and the heading angle tfi, become , 


dft/dt _ _ tan^ 
d^/dt tan^ 


(9a) 


integration of which yields, 


cos^cos^ = cosi 


(9b) 


where, i is the orbital inclination. For small values of cross range angle <p. 
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the orbital inclination i is given by the heading angle itself. Thus, the 
total change in the heading corresponds to the change in orbital inclination 
(plane change). 


Boost and Reorbit 


During the atmospheric flight, the vehicle performs the desired plane 

change and dissipates some energy due to atmospheric drag. Therefore, a second 

impulse is required to boost the vehicle back to orbital altitude. The vehicle 

exits the atmosphere at point F, with a velocity V f and flight path angle 

The additional impulse AV , required at the exit point F for boosting into an 

b 

elliptic orbit with apogee radius R and the reorbit Impulse AV required to 

c c 

insert the vehicle into a circular orbit at point C, are obtained by using the 
principle of conservation of energy and angular momentum at the exit point F, 
and the circularization point C. Thus, we have, 11 

(V +AV ) 2 /2 - p/R = (V -AV ) Z /2 - p/R (10) 

f b ' a c c ' c 


(V +AV )R cosy = R (V -AV ) 

f b a f c c c 


( 11 ) 


Solving for AV and AV from the above equations (10) and (11), 
b c 


AV = 2p( 1/R -1/R )/ri-(R /R ) 2 cos 2 y 1 - V 

b H a c ' L ac f J f 


( 12 ) 


AV = ^ p/R - 2p( 1/R -1/R )/[(R /R ) 2 /cos 2 y -l] 

c ct a c u c a f “ 


(13) 
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Finally, the vehicle is in a circular orbit (of radius R ) moving with 


the velocity = Vp/R^. 


III. OPTIMAL CONTROL 

For the minimum-fuel maneuver, the objective is then to minimize the 
total characteristic velocity for a specified change in heading angle. A 
convenient performance index is the sum of the characteristic velocities for 
deorbit, boost, and reorbit. Thus, 


J - AV + AV + AV 

d b c 


(14) 


Where, AV . AV , and AV sire the deorbit, boost, and reorbit characteristic 

d b c 

velocities respectively, suid su'e obtained from equations (2) and (11) as 


AV = v'mT - (R /R )V cos ( -y ) 

d d ad© © 


(15) 


AV * Vti/R" - (R /R )(V +AV )cosy 

c c a c f b f 


(16) 


Alternatively, AV , AV and AV sure also given by equations (3), (12), 

d b G 

smd (13) respectively. In the normalized form, the performance index becomes, 


J = Av = Av + Av + Av 
d b c 


(17) 


Av = vT7a - (v /a)cos(-y ) 

d ' d e d © 


(18) 


Av « VI /a - f"(v +Av )/a "Icosy 
c ’ c L fb'c J f 

where, 


(19) 


a=R/fc; a=R/fc; Av = AV /Vp/R ; Av = AV MI7R 

d d' a c c' a d d' a c c' a 
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Let us note that for a given circular orbit of radius R , the impulses 

C 

AV and AV are completely determined by the velocity V and the flight path 

be f 

angle y at the the atmospheric exit. The velocity and the flight path 

angle y^ at the entry point are dependent only on the magnitude of the 
deorbit impulse AV . For a specified atmospheric entry (i.e., for a given 

d 

perigee altitude occurring within the atmosphere), we have a fixed value of 
AV and hence fixed values of entry velocity V , and entry flight path angle 
7 ^ as seen from (1) and (2). Therefore, the optimal control problem for the 
minimum fuel consumption is confined only to the segment of the trajectory 
within the atmosphere. Hence, the performance index (14) is more 
appropriately written as 


J = AV + AV (V , y ) + AV (V ,y ) (20) 

d b f f c f f 

Ideally, as seen from equation (12), the minimum value of boost impulse 

AV is zero, when the exit velocity V is made equal to the perigee velocity 
b f 

of the exit-LEO elliptic orbit. Also, the minimum value of reorbit impulse 
AV is obtained when the vehicle exits with zero flight path inclination y . 

c * 

The first step in the optimization procedure using Pontryagin principle 

2 

is to formulate the Hamiltonian as 
H - X bvsiny + X {- A b(l+T) 2 )5v 2 - ■ b S — - ] • 

h v l 1 (b-l + h) 2 J 


A (a bnSvcoso- ♦ b a eosrJ I 

y\ 2 TFTfhT (b _ 1+h) 2 v J 


(b-l+h) vJ 


, ^ fbvcosysini/Tl 
<fi\ (b-l+h) j 
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+ X 




•A 2 b5i?vsln<r _ bvcosycos^tan#") 
(b-l+h) J 


tp\ cosy 


( 21 ) 


where X's are the costates corresponding to the states. The down range angle 
0 does not enter the right hand side of the equation (5) and hence need not 
be considered for the optimization process. 

The optimal control equations for lift and bank angle are given by 


«# .0- ** = 0 

dij ' dar 


leading to 

7) = C w/C 2vX ; tancr = X, A cosy 

LR ' DO v’ VV 

where 

w = j A J + cy cosy 3 2 


( 22 ) 


(23) 


(24) 


The control C L is bounded by the aerodynamic characteristics of the vehicle. 
Thus, for the constrained control, 


ISJ s c l~, or hi 3 c ». 

The costate (adjoint) equations are given by 


(25) 


dA 

h 

ax 

£ 

< 

1 

dx 

“ ah ; 

dx 

A *. 

dX 

dX. 

£ 

dx 

" d* ; 

dx 


ax. 

3v : 


dX 

dtp 


d \ _ ax 

dx ~ 3y 


(26a) 


( 26 b) 
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Boundary Conditions 


The initial and final boundary conditions are given for the normalized 
altitude h as 

h(r=0) = 1.0, h(x=T ) = 1.0 (27a) 

f 

and for the normalized velocity v, and the flight path angle ? as 

(2-v 2 )a 2 - 2a + v 2 cos 2 y * 0 (27b) 

e d d e e 

[2-(v +Av ) 2 "|a 2 - 2a + (v +Av ) 2 cos 2 y = 0 (27c) 

^ f b J c c fb f 

The equations (27b) and (27c) are obtained by eliminating AV from equations 

d 

(1) and (2) and eliminating AV from equations (10) and (11) respectively. 

C 

The remaining multiplier boundary conditions are obtained from the 
transversal ity conditions on the costates. Thus, the optimization procedure, 
requiring the solution of the state equations (8) and the costate equations 
(26) along with the boundary conditions given by equation (27) leads to a 
nonlinear TPBVP, which can only be solved by numerical methods. 

Multiple Shooting Method 

The multiple shooting method is a powerful method for solving nonlinear 
TPBVP. In solving any boundary value problem with the given initial and 
final conditions, we assume additional initial data and integrate forward so 
that the solution satisfies the given final condition as well. This is also 
called a simple shooting method. However, the convergence of the solution is 
highly sensitive to the assumed initial data. The error due to inaccurate 
initial data can be made arbitrarily small by performing the integration over 
sufficiently smaller subdivided panels within the given interval and thereby 
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leading to the multiple shooting method. Thus, the multiple shooting method 
Is a simultaneous application of the simple shooting method at several points 
within the Interval of Integration. Here, the trajectory may be restarted at 
Intermediate points using new guesses. Jacobian matrices are formed for each 
segment. The resulting iteration scheme, based on reducing all 
discontinuities at internal grid points to zero, leads to a system of linear 
algebraic equations. The corresponding OPTSOL code, developed by Deutsche 
Forschungs-und Versuchsanstalt fur Luft-und Raumfahrt (DFVLR) at 
Oberpfaffenhofen, West Germany, was used for solving the present problem. 

IV. NUMERICAL DATA AND RESULTS 

A typical AOTV configuration 12 with L/D of about 1.5 Is shown In Fig. 4. 
The liquid oxygen Is stored In two separate tanks to provide a tapered nose, 
and Inflated chins are used to continue this tapering along the body. A large 
deployable flap is provided to trim the vehicle at low angles of attack for 
maximum L/D performance. A representative set of numerical values used for a 
complete mission from GEO to SSO at an altitude of 556 km Is given below. 3,6 

C =0.1; K = 1.11; m/S = 300 kg/m 2 

DO 

P 8 = 1.225 kg/m 3 ; p = 3.96772xl0 14 m 3 /sec 2 

0 = 1/6900 m" 1 ; R e = 6356.766 km 

H = 120 km; R = 42240.766 km; R = 6912.766 km 

ad c 

Using the above mentioned data, the optimal solution has the following 
entry and exit status. 

Entry status: H = 120 km; V = 10305.58 m/sec 

e e 



Exit Status: H f = 120 km; V = 7462.35 m/sec 

y = 0.1595 deg; <f> - 11.95 deg 

f f 

= 24.1 deg; total flight time * 520 sec 
Characteristic velocities: 

Deorbit characteristic velocity, AV = 1493.32 m/sec 

d 

Boost characteristic velocity, AV = 490.77 m/sec 

b 

Reorbit characteristic velocity, AV * 124.61 m/sec 

C 

Total characteristic velocity AV * 2108.7 m/sec 

The complete mission in terms of the velocity profile is shown in Fig. 

5. Initially, the vehicle is in a circular orbit at GEO moving at a speed V 

d 

= 3064.82 m/sec. A deorbit impulse AV * 1493.32 m/sec is executed to fly the 

d 

vehicle along the GEO-entry elliptic orbit. The elliptic velocity at the 
deorbit point D is V = V - AV * 1571.5 m/sec. At the atmospheric interface 

d d d 

E of altitude H = 120 km, the vehicle attains an orbital velocity V = 

a e 

10305.58 m/sec. During the atmospheric maneuver, the velocity of the vehicle 
is depleted and the exit velocity is V = 7462.35 m/sec. In order to attain 
the desired SSO altitude H * 556 km, a boost impulse AV = 490.77 m/sec is 

c b 

required at the exit F from the atmosphere. Then the elliptic velocity at the 
exit is V = V + AV = 7953. 12 m/sec. The vehicle travels then along the 

f f b 

exit-SSO elliptic path and has a velocity V = 7451.47 m/sec at the reorbit 

C 

point C. In order to insert the vehicle into a circular orbit at this 
altitude H = 556 km, a reorbit Impulse AV = 124.61 m/sec is imparted. The 

c c • 

vehicle is now in a circular orbit at SSO moving with a speed of V = V + 

C C 

AV = 7576.08 m/sec. 

c 
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Fig. 6(a) shows the time history of altitude. The spacecraft enters and 
exits the atmosphere at the altitude of 120 km. The minimum altitude reached 
is 44.72 km. Fig. 6(b) shows a velocity reduction of 2843.23 m/sec. The 
profile of flight path angle with time is shown in Fig. 6(c). The spacecraft 
enters the atmosphere with an inclination of -6.00 degrees and exits with 
0.1595 degrees. The time history of cross range <f> is shown in Fig. 6(d) which 
has a value of 11.95 degrees at the end of the atmospheric maneuver. Fig. 

6(e) shows the variation of heading angle \f), which shows that the atmospheric 
maneuver provides an orbital Inclination of 24.1 degrees. 

The control history is shown in Fig. 7(a). The vehicle enters the 
atmosphere with maximum lift capability and decreases slowly during the 
remaining flight. Fig. 7(b) shows the variation of bank angle during the 
atmospheric flight. Initially the vehicle enters the atmosphere with a bank 
angle of 144.5 degrees to pull the vehicle into the atmosphere but slowly 
drops to about 75 degrees and maintain at a value of 96 degrees for most of 
the remainder of the flight. The lift-to-drag ratio is shown in Fig. 7(c). 

Fig. 8(a) shows the peak heating rate of 402.64 W/sq. cm. As shown in Fig. 
8(b), the peak dynamic pressure is 80.73 KN/sq. m. 

Fig. 9 shows the successive approximations of the altitude H, during the 
course of 0, 15, and 30 iterations in using the multiple shooting method. For 
the sake of clarity only 4 out of 20 Intervals are shown. The initial 
guessed value for the altitude is 120 km at every interval. It can be seen 
how the initially large jumps at the subdivision points are "flattened out" 
with the increase of iterations. 
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V. CONCLUDING REMARKS 


IX LJthls p aper., vie have addressed the problem of minimization of fuel 
consumption during the atmospheric portion of an aeroassisted, orbital 

transfer with plane change^ The complete mission has required three 

. ''■■■' '■ "• ' 

characteristic velocities, a deorbit impulse at HEO, a boost impulse at the 

i '• f,. > ■ ■- ’ - ■ . 

atmospheric exit, and a reorbit impulse at Lfeo. A performance index has been 
formulated as the sum of these three impulses. Application of optimal control 
principles has led to a nonlinear, two-point, boundary value problem which was 
solved by using a multiple shooting algorithm. The strategy for the 
atmospheric portion of the minimum-fuel transfer is to start initially with 
the maximum positive lift in order to recover from the downward plunge, and 
then to fly with a gradually decreasing lift such that the vehicle skips out 
of the atmosphere with a flight path angle near zero degrees. 
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Fig* 1 Space transportation architecture 
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Fig. 3 Coordinate system 
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Fig. 6(c) Time history of flight 
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Fig. 7(c) 
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Fig. 9 Successive approximations for altitude 



FUEL-OPTIMAL TRAJECTORIES FOR AN 
AERO ASSISTED ORBITAL TRANSFER VEHICLE 


Dr. D. S. Naidu 
ODU Research Foundation 
Spacecraft Control Branch, MS 161 
NASA Langley Research Center 
Hampton, VA, 23665-5225 


TABLE OF CONTENTS 


Page 


TABLE OF CONTENTS i 

LIST OF FIGURES ii 

Chapter 

1. INTRODUCTION 4 

2. BASIC EQUATIONS 6 

3. OPTIMAL CONTROL 11 

4. NUMERICAL DATA AND RESULTS 15 

5. CONCLUDING REMARKS 18 

6 . REFERENCES 18 



ABSTRACT: The fuel-optimal problem in noncoplanar orbital transfer employing 
aeroassist technology is addressed. The mission involves the transfer from 
high Earth orbit to low Earth orbit with plane change. The complete maneuver 
consists of a deorbit impulse to inject a vehicle from a circular orbit to 
elliptic orbit to enter the atmosphere, a boost impulse at the exit from the 
atmosphere for the vehicle to attain a desired orbital altitude and finally a 
reorbit impulse to circularize the path of the vehicle. In order to minimize 
the total fuel consumption, a performance index is chosen as the sum of the 
deorbit, boost, and reorbit impulses. For a typical aeroassisted orbital 
transfer vehicle with high lift-to-drag ratio, the simulations are carried out 
using the industry standard POST program. 

NOMENCLATURE 

AOTV: aeroassisted orbital transfer vehicle 

GEO: geosynchronous Earth orbit 

H : altitude 

HEO: high Earth orbit 

J : performance index 

LEO: low Eanth orbit 

OTV: orbital transfer vehicle 

R : distance from Earth center to vehicle center of gravity 

R : radius of Eanth 

E 

SSO: space station orbit 

V : velocity 

y : flight path angle 

AV : characteristic velocity 
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Subscripts 


a : atmospheric boundary 

c : circularization or reorbit at LEO 

d : deorbit at HEO 

e : entry to atmosphere 

f : exit from atmosphere 

1. INTRODUCTION 

The main function of the space transportation system is to deliver 
payloads from Earth to various locations in space. Until now, this function 
has been performed by various rockets, the space shuttle, and expendable upper 
stages using solid or liquid propellants. In particular, considering the 
economic benefits and reusability, an orbital transfer vehicle (OTV) is 
proposed for transporting payloads between low Earth orbit (LEO) and high 
Earth orbit (HEO). The two basic operating modes contemplated for OTV are a 
ground-based OTV which returns to Earth after each mission and a space-based 
OTV which operates out of an orbiting hanger located at the proposed Space 
Station. 

In a typical mission, a space-based OTV, which is initially at the space 
station orbit (SSO), is required to transfer a payload to geosynchronous Earth 
orbit (GEO), pick up another payload, say a faulty satellite, and return to 
rendezvous with the orbiting hanger at SSO for refurbishment and redeployment 
of the payload. The OTV on its return Journey from GEO to SSO needs to 
dissipate some of its orbital energy. This can be accomplished by using an 
entirely propulsive (Hohmann) transfer in space only or a combination of 
propulsive transfer in space and aeroassisted maneuver in the atmosphere. It 
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has been established that a significant fuel savings and hence increased 
payload capabilities can be achieved with propulsive and aeroassisted 
maneuvers instead of all -propulsive maneuvers 1 . This leads to an aeroassisted 
orbital transfer vehicle (AOTV), which on its return leg of the mission, dips 
into the Earth’s atmosphere, utilizes atmospheric drag to reduce the orbital 
velocity and employs lift and bank angle modulations to achieve a desired 
orbital inclination. Basically, the AOTV performs a synergetic maneuver, 
employing a hybrid combination of propulsive maneuver in space and aerodynamic 
maneuver in the atmosphere. 

It is believed that the concept of aeroassisted orbital transfer opens 
new mission opportunities for the space transportation system, especially with 
regard to the establishment of the permanent space station. Fig. 1 shows the 
space transportation architecture relevant to aeroassist technology. The 

2*7 

optimization of fuel is an important aspect of orbital transfer missions 

In this paper, we address the fuel-optimal problem arising in noncoplanar 
orbital transfer employing aeroassist technology. The maneuver involves the 
transfer from HEO to LEO with a plane change and at the same time minimization 
of the fuel consumption. It is known that the change in velocity, also called 
the characteristic velocity, is a convenient measure of the fuel consumption. 
For minimum-fuel maneuver, the objective is then to minimize the total 
characteristic velocity for deorbit, boost, and reorbit (or circularization). 
For a typical AOTV with high L/D capability, the simulations are carried out 
using the industry standard Program to Optimize Simulated Trajectories 
(POST) 8 . Fuel-optimal trajectories are obtained for the given atmospheric 
entry, orbital inclination, and LEO altitude. 
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2. MISSION DESCRIPTION 


The complete mission from HEO to LEO with atmospheric pass is depicted in 
Fig. 2. It is composed of three impulses: first, a deorbit impulse AV at HEO 

d 

to inject the vehicle into a HEO-entry elliptic orbit, second, a boost impulse 
AV at the exit from the atmosphere for the vehicle to attain sufficient 

D 

velocity to travel along an exit-LEO elliptic orbit, and finally, a 
circularizing impulse AV to circularize the path of the vehicle. 

C 

Consider the different phases of deorbit, aeroassist (or atmospheric 
flight),' boost and reorbit (or circularization). 

Deorbit 

Initially, we assume that the spacecraft is in a circular orbit of radius 
R d , well outside the Earth’s atmosphere, moving with a circular velocity = 
Vfi/R. Deorbit is performed by means of an impulse AV , to transfer the 

d d 

vehicle from the circular orbit to elliptic orbit with perigee low enough to 
intersect the dense part of the atmosphere [Fig. 2]. Since the elliptic 
velocity at D is less than the circular velocity at D, the impulse AV is 

d 

executed so as to oppose the circular velocity V . In other words, at point D, 

d 

the velocity required to put the vehicle into elliptic orbit is less than the 
velocity required to maintain it in circular orbit. The deorbit impulse AV 

d 

causes the vehicle to enter the atmosphere of radius R^ with a velocity V^ and 

flight path angle y . It is known that the optimal -energy loss maneuver from 

© 

the circular orbit is simply the Hohmann transfer and the impulse is parallel 
and opposite to the instantaneous velocity vector. 
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Using the principle of conservation of energy and 
deorbit point D, and the atmospheric entry point E, we 


angular momentum at the 
get 9 , 


V 2 /2 -|i/ R = (V -AV ) Z /2 - p/R 
e 7 a dd 7 d 

R V cos(-y ) = R (V -AV ) 

a e a d d d 

from which solving for AV , we get 

d 

AV = ^R" - 2p( 1/R -1/R )/[(R /R ) 2 /cos 2 y -l] 

d d ^ a d 7 L d a 7 e J 


(1) 

( 2 ) 

(3) 


It is easily seen that the minimum value of the deorbit impulse AV 

dm 

obtained at y^ = 0, corresponds to an ideal transfer with the space vehicle 
grazing the atmospheric boundary. To ensure proper atmospheric entry, deorbit 
impulse AV must be higher than the minimum deorbit impulse AV which is 

d dm 

given by 


AV = fjlAC - 2p( 1/R -1/R )/[(R /R ) 2 -l] 

dm dT a d 7L da J 


(4) 


Aeroassist, Boost, and Reorbit 

During the aeroassist (or atmospheric) flight, the vehicle is controlled 
by bank angle with a constant angle of attack to achieve the necessary 
velocity reduction (due to atmospheric drag) and the plane change. Due to the 
loss of energy during the atmospheric flight, a second impulse is required to 
boost the vehicle back to orbital altitude. 

The vehicle exits the atmosphere at point F, with a velocity V and 

flight path angle y . The additional impulse AV , required at the exit point F 

r b 

for boosting into an elliptic orbit with apogee radius R and the reorbit 

C 
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impulse AV required to Insert the vehicle into a circular orbit at point C, 

C 

are obtained by using the principle of conservation of energy and angular 
momentum at the exit point F, and the circularization point C. Thus, we have, 


(V +AV ) Z /2 - |1/R > (V -AV ) Z /2 - p/R 

r b ® c c c 

(5) 

(V +AV )R cosy * R (V -AV ) 

f b a f c c c 

(6) 


Solving for AV and AV from the above equations (5) and (6), 

b c 

AV = 2p( 1/R -1/R )/[l-(R /R ) 2 cos 2 y ] - V (7) 

b 4 a c'*- ac f J f 


AV = «Tm7r" - J2p( 1/R -1/R )/[(R /R ) 2 /cos 2 r -l] (8) 

c c d a c L c a f ■ 

Finally, the vehicle is in a circular orbit of radius R , moving with a 

C 

velocity V = Vp/R . 

c c 

3. OPTIMAL TRAJECTORIES 

For minimum-fuel maneuver, the objective is then to minimize the total 
characteristic velocity. A convenient performance index is the sum of the 
characteristic velocities for deorbit, boost, and reorbit. Thus, 

J = AV + AV + AV (9) 

d b c 

Where, AV , AV , and AV are the deorbit, boost, and reorbit characteristic 

d b c 

velocities respectively, and are obtained from (2) and (6) as 

AV = Vp/R - (R /R )V cos(-y ) (10) 

d d a d e e 
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AV - v^ITT - (R/R ) ( V +AV ) cosy (11) 

c c a c f b f 

Alternatively, AV , AV t and AV are also given by (3), (7), and (8) 

d b c 

respectively. Let us note that for a given circular orbit of radius R , the 

C 

Impulses AV and AV are completely determined by the velocity V and the 

be f 

flight path angle y at the exit. The velocity V and the flight path angle y 

x e e 

at the entry point are dependent only on the magnitude of the deorblt Impulse 
AV . For a given atmospheric entry (l.e. , for a given perigee altitude 

d 

occurring within the atmosphere), we have a fixed value of AV and hence fixed 

d 

values of entry velocity V , and entry flight path angle y as seen from (1) 

e © 

and (2). The optimization problem Is to minimize the total fuel consumption 
(l.e., maximizing the vehicle mass), for a given orbital inlinatlon, and 
apogee altitude at LEO. 

4. NUMERICAL DATA AND RESULTS 

A typical AOTV configuration with high L/D Is shown In Fig. 3. The liquid 
oxygen Is stored In two separate tanks to provide a tapered nose, and inflated 
chins are used to continue this tapering along the body. A large deployable 
flap is provided to trim the vehicle at low angles of attack for maximum L/D 
performance. A representative set of numerical values used for a complete 
mission from GEO to SSO Is given below 10 . 

weight of the vehicle = 112,625 N 
aerodynamic reference area =30.8 sq.m 
aerodynamic reference length = 15.67m 
gravitational constant of Earth = 3.96772xl0 14 m 3 /sec 2 
radius of Earth = 6356.766 km 
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altitude of atmospheric boundary = 120 km 
altitude of GEO ■ 35884 km 
altitude of SSO = 556 km 

Prescribed conditions have the altitude of perigee in the atmosphere as 58 km, 
orbital inclination to be achieved during the amospheric maneuver as 25 
degrees, and the altitude of space station orbit as 556 km. The control is 
achieved by bank angle modulation using a piecewise linear steering option of 
the POST. A constant angle of attack of 10 degrees corresponding to maximum 
lift-to-drag ratio of 1.724 and the 1976 US standard atmosphere, are utilized 
in the simulation. 

Using the above mentioned data, the optimal solution has the 
following characteristic velocities: 

Deorbit characteristic velocity, AV = 1491.04 m/sec 

d 

Boost characteristic velocity, AV = 24.58 m/sec 

b 

Reorbit characteristic velocity, AV ■ 156.36 m/sec 

C 

Total characteristic velocity AV = 1671.98 m/sec 

The complete mission in terms of the velocity profile is shown In Fig. 4. 
Initially, the vehicle is in a circular orbit at GEO moving at a speed V = 

d 

3064.82 m/sec. A deorbit impulse AV = 1491.04 m/sec is executed to put the 

d 

vehicle along the GEO-entry elliptic orbit. A specific Impulse of 456 sec is 
used for all the charact eristic velocities. The elliptic velocity at the 
deorbit point D is V = V - AV = 1573.78 m/sec. At the atmospheric interface 

d d d 

E of altitude H = 120 km, the vehicle attains an orbital velocity V = 
a e 

10305.93 m/sec. During the atmospheric maneuver, the velocity of the vehicle 
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is depleted and the exit velocity is V = 7898.85 m/sec. In order to attain 
the desired SSO altitude H = 556 km, a boost impulse AV = 24.58 m/sec is 

c b 

required at the exit F from the atmosphere. Then the elliptic velocity at the 

exit is V = V + AV = 7923.43 m/sec. The vehicle travels along the exit-SSO 
r f b 

elliptic path and has a velocity V = 7419.74 m/sec at the reorbit point C. In 

C 

order to Insert the vehicle into a circular orbit at this altitude H = 556 

c 

km, a reorbit Impulse AV = 156.36 m/sec is required. The vehicle is now in a 

C 

circular path at SSO moving with a speed of V = V + AV = 7576. 1 m/sec. The 

c c c 

vehicle takes 18,792 seconds for the GEO-entry portion, spends 470 seconds 
during the atmospheric maneuver, and finally requires 1,931 seconds for the 
exit-SSO portion. The total time taken for the complete mission is 5.8869 
hours. 

Fig. 5(a) gives the time history of altitude. The spacecraft enters and 
exits the atmosphere at an altitude of 120 km. The minimum altitude reached is 
52.07 km. The velocity versus time is shown in Fig. 5(b). The vehicle enters 
the atmosphere with a velocity of 10305.93 m/sec and leaves the atmosphere 
with a speed of 7898.85 m/sec, thus giving a velocity reduction of 2407.08 
m/sec. The profile of flight path angle with time is shown in Fig. 5(c). The 
spacecraft enters the atmosphere with an inclination of -5. 169 degrees and 
exits with 1.8596 degrees. Fig. 5(d) depicts the orbital inclination of 25 
degrees achieved during the atmospheric maneuver. Fig. 6(a) shows the 
variation of bank angle during the atmospheric flight. Fig. 6(b) gives the 
peak heating rate as 730 W/sq. cm. As shown in Fig. 6(c), the peak dynamic 
pressure is 30.33 KN/sq.m. According to Fig. 6(d), the maximum g-load is 4.69. 
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6. CONCLUDING REMARKS 


In this paper, we have addressed the problem of minimization of fuel 
consumption during the atmospheric portion of an aeroassisted, orbital 
transfer vehicle with high L/D capability. The complete mission has required 
three characteristic velocities; a deorblt Impulse at HEO, a boost impulse at 
the atmospheric exit, and a reorbit impulse at LEO. A performance index has 
been formulated as the sum of these three impulses. Fuel optimal trajectories 
have been obtained for the vehicle using POST. Future work is concerned with 
obtaining fuel-optimal trajectories with constraints on heating rate and 
acceleration. 
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b) Time history of velocity 
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Fig. 5(c) Time history of flight path angle 
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Fig. 5(d) Time history of inclination 


o 



21 


Fig. 6(a) Time history of bank angle 
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b) Time 
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Fig. 6(d) Time history of sensed acceleration 


